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£^ ' Single field baryogenesis, a scenario for Dirac leptogenesis sourced by a time- 

VO ' 

dependent scalar condensate, is studied on a toy model. We compare the creation 



of the charge asymmetry by the perturbative decay of the condensate with the non- 



Q_i' perturbative decay, a process of particle production commonly known in the context 

CD : 

of inflation as preheating. Neglecting backreaction effects, we find that over a wide 



parametric range perturbative decay and preheating contribute by the same order 



of magnitude to the baryon asymmetry. 



1. INTRODUCTION 

Models for baryogenesis tie together cosmology and particle physics jlj. The discovery 
of small neutrino masses P] and their explanation via lepton-number violating Majorana 
masses and the see-saw mechanism strongly supports the leptogenesis mechanism f2j. An 
alternative to this scenario is to assume pure Dirac mass terms for the neutrinos, that is the 
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absence of Majorana masses, and to induce an asymmetry between the left- and right-handed 
neutrinos, which is subsequently turned into baryon number through sphaleron transitions. 
This idea is referred to as Dirac leptogenesis jj]. 

In many scenarios for baryogenesis the necessary C and CP violation occur simulta- 
neously, induced by a matrix of Yukawa couplings. An exception is electroweak baryo- 
genesis (5I where in first place a CP- violating axial asymmetry is produced which is 
subsequently not erased via equilibration. Parity violation is then contributed in a second 
step due to sphaleron interactions. 

Also Dirac leptogenesis relies on the transition of an axial asymmetry into baryons 
through the sphaleron process. The initial asymmetry is stored within Dirac neutrinos 
and does not get erased due to equilibration until electroweak symmetry breaking since the 
Yukawa coupling to the Standard Model Higss field is tiny. While the CP asymmetry can be 
provided through a matrix of Yukawa couplings and the out-of equilibrium decay of heavy 
scalar particles into Dirac neutrinos it has been suggested that also a single Dirac mass 
term can source CP, provided it is time dependent. This mass term can arise due to Yukawa 
couplings of neutrinos to a rolling scalar field, and the resulting mechanism has been named 



single field baryogenesis p|. Recently, an interesting realisation of this mechanism through 
a decaying Affleck Dine condensate has been proposed |8|. 

Interpreting the scalar condensate oscillating around zero as a large amount of scalar 
quanta at zero momentum, the axial asymmetry can be generated due to the perturbative 
decay of these particles, a 8 connnonly assumed in scenarios for Afnec k -D,ne baryogenesiJl. 
However particles can also be produced nonperturbatively, as first pointed out in Ref. [KJL a 
process which is often referred to as preheating in the context of the decay of the inflaton . 
Various aspects of preheating from the decay of flat directions are discussed in jl^ . 

In parallel, in the coherent baryogenesis [ljj scenario, the oscillating condensate leads 
directly and at tree-level to the production of a charge asymmetry during preheating when 
it couples to matter such that a time- dependent C and CP violating mass matrix arises. 
Consequently, in the case of a single time-dependent mass term a preheating process can 
lead to an axial asymmetry Q] and thereby source single field baryogenesis. 

In the present analysis, we focus on the importance of non-perturbative contributions 
to the baryon asymmetry, and we choose the single field model due to its simplicity. We 
emphasise nonetheless that nonperturbative particle production may be of relevance for 
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other scenarios, e.g Affleck-Dine baryogenesis. 




2. PERTURB ATIVELY SOURCED SINGLE FIELD BARYOGENESIS 

Let us begin by considering a simple toy-model potential 

v = y (l^«l 2 + l^l 2 ) + *y + + X - L( t>^R + KU>R ■ (i) 

The fields <J) U and 4> d are scalar and are multiplets of the electroweak group G E w — SU(2) L x 
U(f)y, 4> u = (2, |), 4>d= (2, — while L = (2, — |) with the components 

/ 

(2) 

and vr = (1, 0) are Weyl fermions. The scalar mass eigenstates are then (3[0 U ] + S[0d|) 
and ^ ( — + both with mass V ^ 2 — m 2 , and eigenstates, (S[<j> u ] + ^[0^]) 

and ^ (9ft[0u] + 3?[0J) with mass V '/i 2 + m 2 . The inflationary Hubble rate is given by H I: 
and we assume that both of these mass eigenvalues are slightly below this value. Therefore, at 
horizon exit the scalar fields get amplified up to a magnitude ~ Hi and a random direction in 
SU(2)i-space. Inflationary expansion then leaves behind a homogeneous vacuum expectation 
value for the scalar fields in our patch of the Universe. This induces large neutrino masses, 
such that they initially do not thermalise. In turn, the potential for <p Ut d does not get altered 
by thermal corrections. 

Since we are interested in neutrino production, in the following, we consider only the 
neutral components 4>° ud . In particular, since the asymmetry is produced from 0°, we take 
this to be the source field as in the single field baryogenesis scenario. Furthermore, we 
assume m <^ ji. 

Coherent oscillations begin at the time when the Hubble rate has decreased to the value 
/i, and the solution for 0° can be approximated for small small /i/H by 



sft[0°] = A? cos ^/i 2 - m 2 ?j + A R cos (^/i 2 + mH 



a-W(t) , (3) 



A[ cos (yV 2 - m 2 ?j + A{ cos (yj 'fi 2 + m 2 ?j cT 3/2 (t) , 



where the values of Af'2 are random initial values arising from inflation, as described above, 
and a(t) denotes the scale factor of the Universe, t denotes comoving time. In order to keep 
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the present discussion simple, we assume Af = A R , = 0, A\ = A 1 and A[ = in our 
patch of the Universe. Under these conditions, the charge density carried by the field (fP u is 
given by 



where we expanded in /i/m and neglected time derivatives acting on the scale factor. Ac- 
cording to the interaction term with the leptons in the potential (0), is transferred to a 
charge asymmetry within the left handed neutrinos when <p° u decays. Due to conservation of 
total lepton number, a precisely opposite amount of the asymmetry is stored within the right 
handed neutrinos. However, this asymmetry in the right-handed sector is not transferred 
into baryons by sphalerons due to the left-handed nature of interactions. 

We assume that the Universe is radiation dominated when coherent oscillations commence 
and that this remains so until the scalar fields decay, such that they contribute only negligibly 
to the entropy density s. Just like Q^, s scales down as a -3 , such that we find for the 
asymmetry within left-handed neutrinos at the time r -1 , when the scalar field decays, 



Here, we have used the relations H = 1.66g* T 2 /mpi , s = ^-g*T 3 and have taken H rs //. 
The number of relativistic degrees of freedom is denoted by g*, such that a is a numerical 
constant of order one for realistic values of g*. Furthermore, we have assumed that 0° decays 
at a total rate T, whereas the decay rate into neutrinos is given by T v = A^/x/(87r), such 
that a branching factor of T u /T arises. 

3. NONPERTURBATIVE SOURCE 

Following Ref. Q|, we calculate the axial asymmetry induced by the nonperturbative 
decay of 0°. We do so by solving numerically the conformally rescaled Dirac equation 







i^ — rriR + i^mj ip — . 



(6) 



We take 




(7) 
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such that by the potential (JIJ, we have 

m R = a\ v ^i[4>% mi = a\ u ^s[(j) u } . (8) 

Furthermore, we introduce the conformal time 77, which is related to comoving time as 
dt = adrj, and we take 3q = d v . 

Let us introduce the positive and negative frequency mode functions, Uh(k, rj) and 
f/i(k, rj) = —i'-f 2 (uh(\s.,i]))*, respectively. They form a basis for the Dirac field, 



^( 3; ) = / 7^4 E e_ik ' X ("^(k) +^l(-k)) , u h =( Lh ) (9) 



where £/, is the helicity two-eigenspinor, h£h = h^. The Dirac equation then decomposes 
into 

id v L h - h\k\L h = m R R h + imjR h , (10) 
i<9,,i4 + h\k\R h = m R L h - imjL h . 

From Lh and Rh, we can define the quantities 

foh = \Lh\ 2 + \Rh\ 2 , hh — \Rh\ 2 — \Lh\ 2 , (11) 
fih = -2n(L h R* h ), f 2h = 2<Z(L* h R h ), 

where f h is the charge density, f 3 h the axial charge density, fih the scalar density and f 2 h 
the pseudoscalar density. Note that one can easily show that foh is conserved by Eq. (jIDj). 
reflecting the charge conservation of the Dirac neutrinos. 

The initial conditions corresponding to a particle number n^k) = \/3o\ 2 are 

/ ')o/. ; ; • ]„/.,, . 

V'k = , |ao| + |Pq| = 1 , (12) 



y ®oRh + PoRh 



where 



oj(k) + hk . m 

i-iu = — 1-, 



\ 2u{k) ' h ~ \m\ \ 2u(k) 



uj(k) — hk 



+ m \m\ 



(13) 



2w(k)(w(k) + hk) ^ 2u(k)(tu(k) - hk) 

and oj(k) = \Jk? + \m\ 2 . Since we assume to have initially zero neutrinos, we take (3q = in 
the following. 
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When (f>1 ceases to oscillate, the particle number is given by 

w&(k) = (hkf 3 h + m R f lh + mifzh) + -. (14) 

Of course there is no charge asymmetry, since there is an opposite amount of antiparticles. 
However, when mj ^ 0, an asymmetry in the number of particles with positive (h = +) 
and negative (h = — ) helicity may be generated. Note that in the limit m^mi — ■> 0, 
nith = \hfzh + §j since then chirality and helicity coincide. Therefore, 

2(n + -n_) = / 3+ + / 3 - (15) 

when the masses vanish. The factor two on the left hand side occurs because the total axial 
asymmetry gets contributions from particles and antiparticles, while /^(k) counts just the 
particles. 

With a prime denoting a derivative w.r.t. rj, the scalar equation of motion reads 

0" + 2-0' + a 2 ^ + ar0 / = O. (16) 

a a<p 

During radiation expansion, a = a R rj, and when H = a! jo? <C \J [i l ± m' 2 , the solution to 
this equation is well approximated by 



A R cos C\J jj? — m 2 -^-rf^j + iA 1 cos (^J fx 2 + m 2 -^rj 2 ^j 



a Sl ))- 3/ Vj rflrf , (17) 



with the same assumptions for the real and imaginary parts as in the previous section. We 
use this solution to obtain the Dirac neutrino mass term (JHJ) and numerically solve Eq. (|10|) 
by integrating up to the time when V > H, such that the Dirac mass term ceases to oscillate 
and the axial charges /3h(k) get frozen in. A typical plot of the spectrum of the generated 
charge charge asymmetry is given in FIG. H Particle production occurs at a time tn es when 
the fermionic mode is in resonance with the coherently oscillating field. The production of 
the soft modes with small momentum k is suppressed because the initial charge asymmetry 
in the scalar field is small {m 2 / fj)t-^ es 1, cf. Eq. (jSJ). Consequently, the production of 
asymmetry within modes which resonate later becomes stronger first. Eventually, there is 
a damping effect due to the red-shifting of the oscillating condensate and finally due to its 
decay at the rate T. Note that due to Pauli blocking —2 < /3+(k) + /3-(k) < 2. 

Of course, the axial asymmetry vanishes in the case when the scalar charge (J3J) is zero. 
When A 1 = 0, the term oc 7 5 in the Dirac equation (JHJ) vanishes and there is obviously no 
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FIG. 1: The axial asymmetry plotted over momentum k, which is taken to be the physical 
momentum at the time when coherent oscillations begin. The choice of parameters is: X u = 0.2, 
A R = A 1 = 20/x, m = 0.05^, T = 0.002/x. 

CP-violation. When A 1 ^ but m = there is also zero scalar charge. Since then the phase 
is constant, d v axg(m r + imj) = 0, the 7 5 -term can in principle be removed at all times by a 
rephasing of the fermionic field. Consequently, even if we do not perform this rephasing, we 
expect to find /3+(k) + /3-(k) = 0, which can also be verified numerically. 

The axial charge density stored within the neutrinos is is the integral over the asymmetry 
within the modes 

7^(/3+(k) + / 3 -(k)), (18) 

where k is to be understood as the physical momentum at the time when coherent oscillations 
begin. The axial asymmetry to entropy ratio then turns out to be 

n{u L ) - n{u L ) _ Q A 

"(/imp,) 3 / 2 ' 1 j 

which we want to compare with the perturbative result (jo]). 

We denote the axial densities n{v L ) — n(u L ) by p res for the nonperturbative or resonant 

case of Eq. (fTHj) and by p pert for the perturbative decay as expressed in Eq. (J3J). The initial 

amplitudes of the scalar field are chosen to be A R = A 1 . We display the produced axial 
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asymmetries over the initial amplitudes in FIGs. E] and EJ where we have taken different 
values for the damping rate T. 

Clearly, the perturbative source p per t gets enhanced by the factor T u /T in (jHJ) as the 
damping T becomes smaller. Note that we have consistently chosen T u < T with the case 
T v = T displayed in FIG. El But also the nonperturbative contribution grows for smaller 
decay rates, because coherent oscillations last longer and a larger phase space volume may 
be filled as the fermionic modes are red-shifted. However, while initially p pcrt and p res grow 
as the square of the initial scalar amplitude, p res gets suppressed for large amplitudes due to 
Pauli blocking, which we do not take into account in our formula for the perturbative asym- 
metry (0). In either case as displayed in FIGs. El and El we note that the total asymmetry 
is the sum of the individual contributions, p per t + Pres- 




5 10 15 20 25 



A n /^=A'/(x 

FIG. 2: The axial asymmetry plotted over the initial amplitude of (j)^ = A R + iA 1 . The choice of 
parameters is X u = 0.2, m = 0.05//, T = 0.01/i. 

4. CONCLUSIONS 

Nonperturbatively sourced single field baryogenesis is a viable scenario. We have shown 
that it contributes over a wide range of parameter space by the same order of magnitude 
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FIG. 3: The axial asymmetry plotted over the initial amplitude of ^ = A R + IA 1 . The choice of 
parameters is X u = 0.2, m = 0.05^, T = T u = 0.0016^. 

as the perturbative source to the baryon asymmetry. Note that a rescaling of the Yukawa 
coupling Aj, can be absorbed into different intitial amplitudes A^ 1 , such that the effect can 
be read of from FIGs. [2]and|21 Besides the model presented here, the coherent baryogenesis 
mechanism is an example for generating the bayon asymmetry directly from preheating. We 
conclude that processes of nonperturbative particle production may be of importance for 
explaining the baryon asymmetry of the Universe. 
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